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Abstract 

We investigate the cosmological evolution of mimetic matter model with arbitrary scalar po¬ 
tential. The cosmological reconstruction , which is the way to construct a model for arbitrary 
evolutions of the scale factor, is explicitly done for different choices of potential. The cases that 
mimetic matter model shows the evolution as Cold Dark Matter(CDM), wCDM model, dark matter 
and dark energy with dynamical Om(z) , where Om(z ) = [( H{z)/Hq ) 2 — 1]/[(1 + z ) 3 — 1], or phan¬ 
tom dark energy with phantom-non-phantom crossing are presented in detail. The cosmological 
perturbations for such evolution are studied in mimetic matter model. For instance, the evolution 
behavior of the matter density contrast which is different from usual one, i.e. 5 + 2H5 — n 2 p5/2 = 0 
is investigated. The possibility of peculiar evolution of 5 in the model under consideration is shown. 
Special attention is paid to the behavior of matter density contrast near to future singularity where 
decay of perturbations may occur much earlier the singularity. 

PACS numbers: 


* E-mail address: jmatsumoto@kpfu.ru 


1 



I. INTRODUCTION 


The existence of dark energy and dark matter is indirectly shown by the observations of 
Cosmic Microwave Background (CMB), Large Scale Structure (LSS) of the Universe, Super 
Novae of Type la (SNIa), and related observational probes. Dark energy is usually regarded 
as a cosmological constant A in the standard model of cosmology, ACDM model. However, 
in general it represents some energy which have a negative pressure with equation of state 
parameter being less than —1/3. On the other hand, dark matter, in particular Cold Dark 
Matter (CDM) means the nonrelativistic unknown matter. There are many candidates for 
dark energy (for review, see 1, 2]) and dark matter 3, J]. However, it is not yet clear which 
model is viable one. 

In this paper, we will consider mimetic matter model j||, which was first proposed as a 
model of dark matter. However, it was realized later that it can be treated as a model of 
dark energy by adding the potential term of the scalar field [§]. Mimetic matter model is 
a conformally invariant theory by making the physical metric as a product of an auxiliary 
metric and the contraction of an auxiliary metric and the kinetic term of the scalar field. 
Whereas, the equivalence between mimetic matter model and a scalar field model with a 
Lagrange multiplier jj] is shown in 8, 9]. Therefore, we treat mimetic matter model as a 
scalar field model with a Lagrange multiplier in this paper. 

The contents of the paper are the following. In Sec. II, we consider the background 
evolution of the LIniverse in general mimetic matter model. The cosmological reconstruc¬ 
tion of the model , which is the way to construct the action which may realize the given 
evolution history of the universe by the choice of potential, is done and ACDM-like era is 
explicitly executed for corresponding scalar potential. The evolution of the matter density 
perturbation is discussed in Sec. III. First, mimetic matter model without potential is used 
to evaluate whether mimetic dark matter also behaves as dark matter at the cosmologi¬ 
cal perturbations level. Next, we investigate the matter density perturbation in general 
mimetic matter model with scalar potential. In Sec. IV, the possibility of the phantom - 
non-phantom transition is explicitly demonstrated in mimetic matter model. The behavior 
of the cosmological perturbations is also discussed, especially the behavior close to finite¬ 
time future singularity. Concluding remarks are given in Sec. V. The units of kp — c — h — 1 
are used and gravitational constant 8ttG is denoted by k 2 = 8tt/Mp\ 2 with the Planck mass 
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of Mpi — G 1 / 2 = 1.2 x 10 19 GeV in this paper. 


II. THE ACCELERATING EVOLUTION OF THE ISOTROPIC AND HOMOGE¬ 
NEOUS UNIVERSE 


We consider the action of mimetic matter model 
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Q: 


S — J d 4 x\/^g 


2k 2 


R-V(<t>) + A(^<R0 + 1) 


5, 


matter ? 


( 2 , 1 ) 


where V is an arbitrary function of the scalar field 0 and A is a Lagrange multiplier. The 
Einstein equation obtained from Eq. (12.1|) is 

Rftv - \diiuR = —K 2 gnvV (</>) - 2k 2 X + K 2 g pu X(g pa d p (j)d a (j) + 1 ) + n 2 T pv , ( 2 . 2 ) 


where R pu = d a T pi/ — d p Y r l (T + T^ and T pu is the energy momentum tensor of 

the usual matter. On the other hand, the equation given by the variation with respect to 
the scalar held </> is the following one: 


- V* - 2W i {Xd li <j>) = 0, (2.3) 

where = dV((j))/d(j). The constraint equation given by the variation of A is 

g^d^dvfy + 1 = 0. (2.4) 

When the spatially hat Friedmann-Lemaitre-Robertson-Walker (FLRW) metric, ds 2 = 
— dt 2 + a 2 (t) Y^=i dx l dx l ), is taken, FLRW equations are written by 

3 H 2 = k 2 P - 2k 2 A0 2 + k 2 V + k 2 A(</> 2 - 1), (2.5) 

-2 H - 3 H 2 = k 2 w P - n 2 V - k 2 A(0 2 - 1), (2.6) 

where H = d[t)/a{t). p is the energy-density of the matter and w is the equation of state 
(EoS) parameter expressed by w = p/p. From Eqs. (j2.3(1 and (12.41) we have 

2A0 + 2A0 + 6LTA0 — = 0, (2.7) 

0 2 = 1 . ( 2 . 8 ) 

The equation of continuity of the matter is 

p + 3(1 + w)Hp = 0. (2.9) 
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If w is constant, one can integrate Eq. (12.91) as 


P = Po\ — 
a 0 


—3(l+ui) 


where poa^ 1+w ^ is an integration constant. While, combining Eqs. (j2.5|) . 
gives the equations 


and 


A — —-zH + —(1 + w)p 


K 


and 


V = wp+ —(2H + 3H 2 


K 


( 2 . 10 ) 


it 


( 2 . 11 ) 


( 2 . 12 ) 


Therefore, all the solutions for the variables will be found if Eq. (12.121) can be solved with 
respect to a(t). From the other side, one can say that this model is able to realize arbitrary 
evolution history of the Universe by using the arbitrariness of the function V{(j>). Indeed, 
one can choose the form of V(<p) for the arbitrary function / such that a(t) = f(t ) as 

— 3 ( 1 + 1 /;) 

- 4 - - 

-2 


vw=Wf, ° [ .m) 


1 ( J'\4>) , fW 

« 2 V m /(0) 2 


(2.13) 


A. Cosmological reconstruction 

Some forms of the function V(4>) are considered in ref. 6]. Here, we consider some other 
examples of V (</>), which are adequate to the cosmological observations and yield the acceler¬ 
ating universe. In particularly, the form of V (0) which can yield the phantom - non-phantom 
transition is treated in Sec. IIV1 In the following, we assume the value of EoS parameter of 
the usual matter as Pmatter/Pmatter = ^matter = 0. 


1. Mimetic wCDM model 

wCDM model consists of cold dark matter and a fluid with EoS parameter w < 0. For the 
case w = —1, wCDM model is equivalent to the ACDM model, so its background evolution 
can be realized by choosing the form of the potential as V = A. In the case of w ^ —1, 
we can reproduce the backgroimd evolution of the Universe in wCDM model if the effective 
energy density —2A + V and the effective pressure — V satisfy the following conditions: 

- 2A + V = Pdmoa- 3 + p w0 a~ 3{1+w) (2.14) 
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and 


* V = wp w0 a~ 3(1+w \ (2.15) 

where pd m o and Pw o are the current energy-density of dark matter and that of a fluid, 
respectively. Here, the current value of the scale factor is chosen as a 0 = 1. Thus, 

A = { Pdmoa ~ 3 + (1 + w)p w0 a~ 3(1+w) j , (2.16) 

H = -wp w0 a~ 3(1+w) . (2.17) 


The specific form of the potential is very complicated even if it can be obtained by solving 
Eq. (12.51) . However, one can easily obtain some approximated forms of V{(j>). 

When w > — 1, an approximated solution of wCDM model is given by 


a(t) = 


^ ^ ( Phi) T pdm {)) 


2 
1 3 


1 + 


(1 + wf 


PwO 


PbO + PdmO 


. K {.PbO T Pdm,o)t 


3(1+™) 


(2.18) 

where pbo is the current energy-density of the baryonic matter, because the solutions of 
Eq. (12.51) in the matter dominant universe and in the fluid matter dominant universe are given 
by a = (3fi; 2 (pbo+ Pdmo)/4) 1 ^ 3 t 2 ^ 3 and a = [3(1 + w) 2 K 2 p w0 /^\ 1 ^ 3( ' 1+w ' > h 2 ^ 3( - 1+w ' ) \ respectively. 
Therefore, an approximated form of potential V{4>) is expressed as 


v(4>) = 


9k 2 (J) 2 


1 + 


(1 + w ) 2 


Pw 0 


PbO+pdmO 


(|« 2 (PbO + PdmO ) 0^ ) 


'3(1+™) 


X 


3 w + w 2 


+ 


w 


(1 + w) 2 (1 T w') 2 


1 + 


3(1+™) 


. (2.19) 


f 1 oW)' 

Eq. (I2.18P is a very rough approximation as a kind of the interpolation. To improve the 
accuracy of the approximation, one can use the current value of the Hubble rate H 0 as: 

"3 


a(t) = 


" (PbO T f plrn li ) t 


+ 


-(1 + W) K p w0 t 


3(1+™) 


2-\-w 


+ c(H 0 t)^y, 


( 2 . 20 ) 


where C is a constant tuned to satisfy a/a = Hq when a(t) = 1. The value of C is evaluated 
by solving a simultaneous equation (a(t 0 ) = 1, H(t 0 ) = H 0 } with respect to t 0 and C. In the 
case of w = —0.7, to = 116.752 x 10 8 years and C = 0.126265 are obtained. The comparisons 
between the exact wCDM model and Eq. (12.20ji when w = —0.7 and w = —0.9 are shown in 
Fig. □ It seems that Eq. (I2.20p is not a bad approximation when w is large, although only 
three points t = +0,t 0) +oo are used as the interpolation points. 
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Figure 1: Comparison of the time dependence of the scale factor between wCDM model and 
Eq. (12.201) . Thick solid line, solid line, dashed line, and dotted line represent Eq. (12.201) with 
w = —0.7, Eq. (12.201) with w = —0.9, wCDM model with w = —0.7, and wCDM model with 
w = —0.9, respectively. Hq = 74 km/s/Mpc, fl m o = 0.3, and fi w o = 0.7 are assumed. 

2. Reconstruction with respect to Om 


One can rewrite Eqs. (12. 11 j) and (12.12j) by using Om(z ) and Om'(z) instead of H{t) and 
H'{t) as 


Hi 


K 


A = < -3/2(1 + z) 3 Om - -(1 + z) A Om' + -(1 + z)Om' \ + -p 60 (l + ^) 3 , 


K 2 


2 y 2 

{3(1 — Om) — (1 + zYOm' + (1 + , 

where pb 0 is current energy-density of the baryons. Om(z ) is defined by HQ 

W,)/j/.)»- i 


( 2 . 21 ) 

( 2 . 22 ) 

(2.23) 


(l + ^) 3 -l ' 

In the case of the ACDM model, Om(z ) is in accord with the current matter fraction f2 m o- 


However, in the recent paper[15], a discrepancy between CMB observations (Planck+WP) 
and BAO observations is pointed out. It is pointed out that the value of Omh 2 is larger 
for larger redshift. The explicit values of fl m oh 2 for Planck+WP and Omh 2 for BAO ob¬ 
servations are O m o^ 2 = 0.1426 ± 0.0025 16] and Omh 2 « 0.122 ± 0.01 15], respectively. Of 
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course, these values are incompatible in the ACDM model, because Omh 2 is constant with 
respect to redshift. However, we can construct a model that satishes these conditions. If we 
consider the following scale factor as a solution of Eqs. (j2.5jl and (12.61) : 

1/3 


a(t) = 


9 if) 


1 - 9(t) 
g(t) = (0.122 + 0.04 


smh 2/3 
tcMB 


;\/l - g{t)H 0 t 


0.02^ mb 


/h\ 


(2.24) 


t t 2 

where tcMB is a time of * = 1089, then Omh 2 rs 0.122 for small redshift and Omh 2 = 0.142 
at z — 1089 is realized. The potential which causes such an expansion history is given by 
using Eq. (j2.13j) . however, it is not written down explicitly, being very complicated. Similarly, 
we can reconstruct any requested universe history. In the next section, the evolution of the 
matter density perturbations is studied. 


III. COSMOLOGICAL PERTURBATIONS WITH ACCOUNT OF MATTER 


In this section, we investigate the growth rate of the matter density perturbations (for 


general review of cosmological perturbations theory, see 12 


- 3 ). 


When we use the metric of 


the Newtonian gauge, ds 2 = (—1 + 2&)dt 2 + (l + 2T)a 2 (f) ^ 3 =1 dx l dx\ then (00), (0*) = (z0), 
(ij), i 7 ^ j, and (ii) elements of the linearized Einstein equation in the Fourier space are 
represented as follows, respectively: 


k 2 

6 H 2 $ - 2—T - 6 Hd^ = - k 2 6 P + 2k 2 5\ + 2k 2 \<& + 2k 2 A<50 - n 2 VJ(j), (3.1) 

GT 

2 (H$ + c) 0 T) = /t 2 (p + p)5u + 2k 2 X 5<j>, (3.2) 

$ - + = 0, (3.3) 


k 2 m 
+ ' 


- 2Hdn — AH — 6H 2 ) <h — 


k 2 d l d l 

a 2 a 2 


+ 2d 0 8 0 + 6Hd 0 j T = 


n 2 c 2 6p — 2firA(<l> + 5<j)) — k 2 U 0 < 50 , 


(3.4) 


where we defined the sound speed by c 2 = 5p/5p, the matter density perturbation by 
S = 5p/p and the scalar perturbation of the four velocity as di5u. The linearized equations 
of (12.3p and (I2.4[) are given as follows: 


5X + 3HSX + A Sip + (A + 3H\)6j> + 




50 
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+A$ + 2(A + 3H\)Q + 3AT = 0, (3.5) 

<j)8(j) + 0 2< f> = 0. (3.6) 

From the perturbation of the equation of continuity, = 0, one gets 

<5/5 + 3 H(5p + 5p ) + cT 2 <9j{(p + p)<9j<5n} + 3'F(p + p) = 0 , (3.7) 

GT 3 <9 0 {a 3 (p + p)di5u} + c 2 (9i5p - (p + p)<9i<L = 0 . (3.8) 


In the following, we treat the equation of state parameter and the sound speed of the matter 
as w = c s = 0 by assuming the era after matter dominance. 

A. Mimetic dark matter 

First, we consider the case E'(0) = 0 as a simple example. Then, the field A behaves 
as dark matter in the flat FLRW space-time jsj]. Let us investigate whether the field A also 
behaves as dark matter when we consider cosmological perturbations. First, Eqs. (13. 6 p and 
(13.3p give the relation <3> = T = —<50. By applying the subhorizon approximation, <9 0 ~ H 
and H 2 k 2 /a 2 , Eqs. (13.11) and (13.51) . then, yield the following equations: 

k 2 ■ 

2—<50 ~ — n 2 8p + 2k 2 <5A, (3.9) 

a z 

k 2 

<5A + 3H5\ + A—<50 ~ 0, (3.10) 

a z 

Differentiating Eq. (I3.10p with respect to t gives 

+ 2 H5 X + ^<50 ~ 0, (3.11) 

a z 

where <5 A = 5 A/A. In the case of w = c 2 = 0, we can obtain the following equation by using 
Eqs. (J33D, (USD , 

. k 2 ■ 

<5 + 2H8 -I—-<50 ~ 0. (3.12) 

a z 

Taking into account p, A oc a -3 , we obtain 

Stot + 2H5 tot + ^<50 ~ 0, (3.13) 

a z 

where p tot = p — 2A and <5 to t = (8p — 2<5A)/(p — 2A). By using Eq. (13.9p . we finally obtain 










Now, one can recognize the equivalence between mimetic dark matter and fluid dark matter 
at the linear perturbation level because Eq. (13.14p is same as that of the ACDM model. The 
equivalence of mimetic dark matter and fluid dark matter had been shown in Ref. [7]. 

B. Mimetic matter with potential 


In this section, we consider more involved case V^ 7 / 0. The perturbation equations (13.21) . 
(13.4j) . and (13 . 7|) can be simplified as follows by using Eqs. (13.4[) and (13.6[) : 


—2(50 — 27/(50 = n 2 p5u + 2k 2 A (50, 

2(50 + 8//<50 + (477 + 6// 2 )<50 = -kV 0 (50, 

k 2 

Sp + 3H5p -— 3pScj) = 0 . 


(3.15) 

(3.16) 

(3.17) 


Eliminating the term proportional to 5u from Eqs. (13.15ft and (13.17ft gives 


r . 2 k 2 \ ... 277 k 2 ■ 2A fc 2 n 

5 + ( —-—- — 3 ) (50 + ———^0 + ~^y^0 — 0. 


K 2 p a 2 


n 2 p a 2 


(3.18) 


Then, we can derive the following evolution equation of the matter density perturbation: 




1677 2 - 4 K 2 (p - 2A) - 577- + O 


a 2 H 2 


- ( 877 3 - 3n 2 H(p - 4A) - 477 2 ^ + n 2 p'~ + O 
2 \ A A 


A ' “ V k 2 
A „ / a 2 H 2 


k 2 


5 = 0. 


(3.19) 


Eq. (13.191) is quite different from that in quintessence model. There are no oscillating solu¬ 
tions, moreover, the quasi-static solutions are modified. In fact, if we assume 

-2 


K 


then 


<5 + 2/7(5 - — p6 = 0, 


777 - - j 5 + { 1677 2 - 4 n 2 {p - 2A) - 5Hj ) 5 

4 f 877 3 - 3k 2 H( P - 4A) - 477 2 ^ + k 2 P ^ ) 6 - ^-p\6 = 0. 
2 \ X A / 2 


(3.20) 


(3.21) 


Therefore, the solution which satisfies Eq. (13.20)1 is not the solution of Eq. (13.19)1 except for 
the cases p = 0 or A = 0 or <5 = 0. It is necessary to be careful because there is an ambiguity 
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to express Eq. (13.2ip due to adding the terms proportional to Eq. (13.20P or the derivatives 
of Eq. f!3.20p are allowed. In /c-essence model, the growth equation of the matter density 
perturbation is also represented by fourth order equation[l7] as in Eq. (13.191) . However, in 
general, the sound speed read off from the field equation (13.51) is not zero in the fc-essence 
model, so there are no modifications to the quasi-static solution. Thus, the modification to 
the quasi-static solution in Eq. (13.19P is the peculiar feature in this theory. However, in the 
matter dominant era or in A dominant era, the term —K 4 pA5/2 is less than the other terms 
in Eq. (13.2ip . Therefore, we need to consider the era where there is almost same amount 
of dark energy and dark matter such as the current universe to distinguish Eq. (13.19p from 
Eq. (13. 21 p . Whereas, in the case of V0, = 0, Eq. (I3.19p recovers the differentiated equation 
of (I3.14p . correctly. 

The procedure to obtain Eq. (13.19P is 


1. differentiating Eq. (I3.18p with respect to t 

2. eliminating the terms proportional to 5 0 from the equation derived at the first step 
by using Eq. (13. lfij) 

3. differentiating the equation derived at the second step with respect to t 

4. eliminating the terms proportional to 5 0 from the equation derived at the third step 
by using Eq. (13.16P 

5. differentiating the equation derived at the fourth step with respect to t 

6. eliminating the terms proportional to 50 from the equation derived at the fifth step 
by using Eq. (I3.16P 

7. eliminating the terms proportional to 50 from the equations derived at the second, 
fourth, and sixth step by using Eq. (I3.18j) 

8. eliminating the terms proportional to 50 and 50 by using the equations derived at the 
seventh step. 
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Figure 2: Redshift dependence of the growth rate function f(z) in mimetic wCDM. The dot-dash 
line, the solid line, and the dashed line express f(z) in the cases w = —1.2, w = — 1 , and w = — 0 . 8 , 
respectively, = 0.73, f ^ m0 = 0.23, and fl&o = 0.04 are assumed to depict the above figure. 
Initial conditions are f"( 100) = /'(100) = 0 and /(100) = 1. 

1. Mimetic wCDM 

Let us evaluate the difference between Eqs. (13.19jl and (13.20)1 in the case that mimetic 
matter behaves like wCDM model. For a correct evaluation, we will not use Eq. (12.1911 but 
use Eqs. (12.161) and (12.171) . If we utilize N = In a instead of time t, we can describe all 
the functions in Eq. (13. 19ft as the functions of a by using Eqs. (12.51) . (12.61) . (12.lOj) . (12.161) . 
and (12.171) . Therefore, Eq. (13.191) can be solved at least by numerical calculations. Fig. [2] 
shows the changes of the growth rate function f(z) = dlnS(z)/dN, with respect to z, where 
z = —1 + 1/a = —l + e~ N . Figs. [3] and 0] express the differences between mimetic wCDM and 
wCDM model in the growth rate function. Here, = 0.73, O dm0 = 0.23, and Qbo — 0.04, 
where = k 2 pxo/{3Hq) is the current energy fraction of the matter X, are assumed to 
depict the figures. As the initial conditions, /"( 100 ) = /'(100) = 0 and /(100) = 1 are used, 
because, in the matter dominant era, i.e. V — const., we could use Eq. (13.141) and can obtain 
f"(z) = f(z ) = 0 and f(z ) = 1. As noted in Sec. IIII Al and at the beginning of Sec. IIIIB1 
there is no difference between mimetic model and wCDM model when w = — 1. However, 
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Figure 3: Difference between mimetic wCDM and wCDM model in the growth rate function f(z) 
when w = —0.8. 


in the case w ^ — 1 & w < 0, there is difference which come from the term —K 4 pXS/2 in 
Eq. (13.211) . One can see in Figs. [3] and |4] that the differences are small but surely exist. The 
term —K 4 pX5/2 is not relevant at the matter dominant era, so we can also find that the 
differences become small as redshift. z becomes large. 

2. Dynamical Om(z) 

The behavior of the growth rate function f(z) when the background solution of the 
Universe is given by Eq. (12.241) is expressed in Fig. [5] Hq = 100 hkms~ 4 Mpc^ 1 , h = 0.7 is 
assumed to depict Fig. [5] Although Om(z) is a smoothly changing function from 0.142/h 2 
to 0.122/h 2 when we consider z < 1089, however, the growth rate function f(z) is larger 
than that of the ACDM model with Q m0 h 2 = 0.142. This effect comes from the difference 
between Eqs. (I3.19P and (13.201) . In fact, if we evaluate f(z) for a constant Om model with 
Omh 2 = 0.142 in Eq. (13. 19|) . we have a larger value than that of the dynamical Om(z ) 
CUED. 
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Figure 4: Difference between mimetic wCDM and wCDM model in the growth rate function f{z) 
when w = —1.2. 



z 


Figure 5: The growth rate function f(z) when the evolution of the background space-time is 
described by Eq. (12.241) . Solid, Dashed, and Dot-dashed curves represent f(z) in this model, f(z ) 
in the ACDM model with D m o/i 2 = 0.142, and f(z ) in the ACDM model with D m o^ 2 = 0.122, 
respectively. Initial conditions are /"(100) = /'(100) = 0 and /(100) = 1. 
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IV. ON THE POSSIBILITY OF THE PHANTOM - NON-PHANTOM TRANSI¬ 
TION 


A. Reconstruction of the background evolution 


In the quintessence model, the value of EoS parameter w = p#/p$ is constrained as 
— 1 < w < 1 so that phantom behavior w < — 1 cannot be realized. Here, we formally 

consider the possibility to construct a model, which can describe the transition from — 1 < w 

to w < —1 or from w < —1 to —1 < w. If the Hubble rate H(t) is given by the following 
function, 

cv cv 

H(t) = - (T 0 + t ) 3 - /3(T 0 + t) + 7 , 7 = ——Tq + /3T 0 , (4.1) 

H(t) = ho + h\ sin(z/f), (4.2) 

(4.3) 


H(t) — c 0 ( - + 


or 


H(t) = t 2 


t t s — t 


+ 


t 2 o~t 2 


t 2 + t 2 


(4.4) 


then for the case, 0 < t < t_ or t + < t, h 1 z/cos(z/t) > 0 , t s /2 < t < t s , or t'_ < t < 0 
or t > t' + , where t± = —T 0 ± y/P/a and t± = (tjj — t\)/ 2 , effective EoS parameter 
w e ff = —2H/(3H 2 ) — 1 satisfies w e g < —1, respectively|18]. Here, T 0 , hi, and u are real 
valued constants and a, /3, h 0 , Co, t s , to, t\ < 1 0 and t 2 are positive constants. Eq. (14.111 could 
also describe the unified model of the inflation and dark energy because there are two regimes 
0 < t < t- and t + < t satisfying w e g < — 1 , the era 0 < t < t_ and the era t + < t correspond 
to the inflation and a dominance of dark energy, respectively. The forms of the potential 
which cause the above expansion rates of the Universe are obtained by using Eq. (12.13(1 as 


m = 


fir 


a 


y (T 0 + 0 ) 6 - 2otf3(To + 0 ) 4 + 2 « 7 (T 0 + 0 ) 3 


+ (2a + 3/3 2 )(T 0 + 0 ) 2 - 6/37(T 0 + </>) + 3^ 2 - 2/3 


V{4>) — — [3 hi + 6h 0 hi sin (up) + 2h, 1 v cos(up) + 3 h 2 sin 2 (zz</>)] , 


V(P) = 


«r 


3c 2 


2cn 


+ 


6 c o 


+ 


2cq + 3c5 


P(t s - 4>) (t s - 4>) 


and 


4 t 2 p + 3 1 2 


+ 


6 1 2 


+ 


3 1 2 


4 t 2 p 


' (<w l (<§ - ifi - (tf + 00 1 <(? + ry j ’ 


(4.5) 

(4.6) 

(4.7) 

(4.8) 
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respectively. Here, EoS parameter of the matter is assumed to be 0. As seen in Eq. (j2. 16jl . 
the phantom crossings in these examples are realized by the transition from A > 0 to A < 0 
when there is no usual matter. Hence, in mimetic matter model with potential the phantom- 
divide transitions maybe realized by proper choice of the scalar potential. 


B. Growth of the matter density contrast in phantom crossing model 


Let us consider the evolution behavior of the matter density contrast in the model with 
Eq. (14.7p as an example. In this model, the Universe will experience the singular evolution in 
H{t ) if t s is larger than the current time, ft is known that there may occur several finite-time 


future singularity scenarios of the Universe as 
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2l| which are classified as [201 ]: 


• Type I (’’Big Rip”) : For t —* t s , a —> oo, p eS —> oo, and |p e ff| —> oo. This also includes 
the case of p e g, p e s being finite at t s . 

• Type If (’’Sudden”) : For t —> t s , a —>■ a s , p e g —> p s , and |p e ff| oo. 

• Type Ill : For t —* t s , a —> a s , p e R —> oo, and |p e frI —>■ oo. 

• Type IV : For t —> t s , a —> a s , p e g —> 0, |p e ff| —>• 0, and higher derivatives of H diverge. 
This also includes the case in which p e ff (Peff) or both of p e g and p e fr tend to some 
hnite values, while higher derivatives of H diverge. 


Singularity in the case under consideration corresponds to Type I or Type III. When c 0 = 2/3 
and t s — 4 x 10 10 years, we have a solution of Eq. (14. 3 K as 

2 

(yb) ’' (4 ' 9) 

where C\ is an integration constant. It is chosen as C\ = 1.6 to comply with the observed 
background evolution of the Universe. In this case, lUgg is simply expressed as ?u e ff = —2 t/t s . 
The evolution behavior of the matter density perturbation is changed around the singularity 
because the approximation \a 2 H 2 /k 2 \ 1 no longer holds in Eq. (13.191b In Fig. [HI 

the behavior of the growth rate function f(z ) = d\\\5(z)/dN(z) until the current time is 
expressed. On the other hand, the future evolution of f{t ) is expressed in Fig.|71 Fig.Hshows 
singular behavior in f(t) around t = 2.7 x 10 10 years, although the singularity in a(t ) and 
H(t ) occurs at t = 4 x 10 10 years. The reason why such a behavior of fit) occurs is instability 
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z 


Figure 6: The growth rate function f(z) in the model (14.71) 


of the perturbations of phantom matter even at the classical level. To be concrete, the flip 
of the sign of the kinetic term would cause the instability. A decay of perturbations much 
earlier of future singularity found in js2| is just the same as this instability (for discussion 


of problems of phantom-divide transitions in perturbations, see 22-24J). 


V. CONCLUSIONS 

In this paper, we have investigated the evolution of the matter density perturbation 
in mimetic matter model. In Sec. II, the reconstruction, the way to construct a model 
which can describe arbitrary evolution history of the Universe, has been explicitly shown 
and some examples, including mimetic wCDM model and evolving Om(z) model, have been 
given. At the beginning of Sec. Ill, the growth rate of the matter density perturbation has 
been investigated when potential term is not dynamical. As a result, it has been shown 
that mimetic dark matter also behaves as CDM at the perturbation level. In the case of 
dynamical potential, it has been shown that the differential equation of the matter density 
perturbation is changed. In general, this equation is different from that in the ACDM model 
and the quasi-static equation of the matter density perturbation in quintessence model. We 
have presented the behavior of the growth rate function in mimetic wCDM model and in 
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Figure 7: The future evolution of the growth rate function f(t) in the model M.7I) 


a model of dynamical Om{z ) as examples. The difference seems not to happen in mimetic 
wCDM, but it clearly appears in dynamical Om(z). In Sec. IV, we have proved the possibility 
to realize phantom — non-phantom transition in mimetic matter model. It has been shown 
that phantom - non-phantom transition is possible because the sign of kinetic term of the 
scalar field is determined by the sign of the Lagrange multiplier. However, usually phantom 
brings an instability. Such an instability occurs not only as quantum effect but also as 


classical effect 
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, and it breaks the perturbation theory as is explicitly demonstrated. 


It would be interesting to study the accelerating universe evolution and cosmological 
perturbations in the generalizations of mimetic matter model which include mimetic F(R) 


gravity 
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36], higher-derivative mimetic model |34J or its further generalizations 35]. In 


this case, we have not only scalar potential but also extra function like F(R) or higher- 
derivative term. Hence, the cosmological reconstruction can be executed quite easily. This 
will be studied elsewhere. 
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